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SUMMARY 


An analytical study was undertaken in order to further 
the understanding of the interaction of airborne and struc- 
tureborne noise radiated by aircraft materials. The theory 
and the results of several computer simulations of the noise 
radiated by thin, isotropic, rectangular aluminum plates due 
to fully coherent combined acoustic and vibrational inputs 
is presented. The most significant finding of the study was 
the extremely large influence that the relative phase 
between the inputs has on the combined noise radiation of 
the plates. It is shown that phase dependent effects 
manifest themselves as cross terms in both the dynamic and 
acoustic portions of the analysis. The computer simulations 
show that these cross terms can radically alter the combined 
sound power radiated by plates constructed of aircraft-type 
materials. The results of the study suggest that airborne- 
structureborne interactive effects could be responsible for 
a significant portion of the overall noise radiated by 
aircraft- type structures in the low frequency regime. This 
implies that previous analytical and experimental studies 
may have neglected an important physical phenomenon in their 
analyses of the interior noise of propeller driven aircraft. 
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1. INTRODUCTION 


High interior noise levels in propeller driven aircraft 
have historically been a cause for concern in both the com- 
mercial and military sectors of the aircraft industry. 

Future problems with interior noise levels in aircraft are 
expected to intensify due to the advanced turboprop prop- 
ulsion systems now being incorporated into the design of 
transport aircraft. 

The noise entering the cabin of a propeller driven 
aircraft is generally divided into two major categories, 
viz. airborne noise and structureborne noise. Airborne 
noise is generated aerodynamically by the propellers, pro- 
pagates along an acoustic path, and is transmitted through 
the sidewalls of the aircraft into the interior. Structure- 
borne noise has its source in the vibration of the wings or 
other structural members of the aircraft. This vibrational 
energy propagates along structural paths into the aircraft 
where it causes vibration of the various surfaces in the 
cabin and ultimately radiates noise. 

The distinction between airborne and structureborne 
noise is important because the methods typically used in 
reducing airborne and structureborne noise are quite dif- 
ferent. For example, if the predominant source of the noise 
in the aircraft is structureborne, then the problem might be 
solved through the use of vibration isolators, or by the 
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application of damping materials. If the predominant source 
of the noise is airborne, then the problem might be solved 
by aft mounting the propellers or by adding massive mater- 
ials to the sidewalls. 

Because of the different treatment methods for the 
airborne and structureborne components of the noise, re- 
searchers, in the past, have devoted their attention almost 
exclusively to either the airborne path or the structure- 
borne path. (See reference [1] for a summary of past 
research on the airborne path and references [2]- [4] for 
some of the more recent efforts aimed at structureborne 
noise.) While this separatist approach may produce more 
efficient methods for treatment of the two individual paths, 
it neglects the possible interaction that may take place 
between the airborne and structureborne components. An 
interactive component arises due to the coherent nature of 
the noise propagated along these two paths. (The airborne 
and structureborne noises are fully coherent since they both 
have the same source, viz. the propellers.) Presently, the 
relative importance of these interactive effects in pro- 
peller driven aircraft is unknown. The purpose of this 
paper is to present the theoretical developments of a study 
[5] which was aimed at furthering the understanding of the 
interaction between the airborne and structureborne noise. 

In addition to the theory, the results of several computer 
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simulations are presented for the noise radiated by aircraft 
type materials due to fully coherent combined acoustic and 
vibrational inputs in the low frequency regime. 

2. PROBLEM APPROACH 

In order to fully understand the interaction between 
airborne and structureborne noise, both analytical and 
experimental studies were performed. The analytical model 
was used to simulate the noise radiation of a structure due 
to simultaneously combined coherent acoustic and vibrational 
inputs. The experimental study was then used to verify the 
behavior predicted by the analytical model. This paper 
presents the details of the analytical model and the results 
of the computer simulations. The apparatus used for the 
experimental study and the results obtained from the experi- 
ments will be presented in another paper. 

Simple isotropic rectangular plates served as the test 
vehicle for the studies. Plates were chosen because they 
possess most of the vibrational and sound radiative pro- 
perties that are exhibited by actual aircraft sidewalls (due 
to their thin shell construction). The analytical model 
assumes that the plates are rectangular and simply supported 
in an infinite, rigid baffle. The sound radiated by the 
plates is assumed to be generated by the flexural (bending) 
vibrational response of the plate and is radiated to a free- 
field acoustic space (anechoic). Futhermore, the dynamic 
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response of the plate (which generates the sound) is assumed 
to be dependent only on the incident airborne and structure- 
borne forcing functions on the plate and independent of 
(uncoupled from) the sound pressure radiated by the plate. 
The plate geometry is defined in figure (1). 

A normally incident, spatially uniform pressure field 
was used to model the airborne input. A point load was used 
to model the structureborne input. The analysis was confin- 
ed to the 0-1000 Hz frequency range since this range encom- 
passes the most troublesome noise region for propeller dri- 
ven aircraft. 

The effects of several parameters on the interaction 
between the airborne and structureborne components were 
investigated. Parameters studied included the relative 
magnitude and phase of the acoustic and vibrational inputs, 
the location of the structureborne input, and the level of 
structural damping. 

3. ANALYSIS OF THE DYNAMIC RESPONSE 
3.1 FREE VIBRATIONAL RESPONSE FOR THE UNDAMPED CASE 

The governing differential equation for the free 
dynamic response of an undamped, simply supported, rect- 
angular orthotropic plate (from reference [6]) is given by 

D n3 4 u/9aJ + 2(D 12 + 2D 66 )3 4 u/3a^3a2 + D 22 3 4 u/3a 4 

+phu=0, (1) 
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where the plate rigidity constants that relate the internal 

f 

bending and twisting moments of the plate to the twists and 


curvatures they induce are given by 

D n = h 3 /12 E 1 /(l-y 1 y 2 ) , (2) 
D 22 = h 3 /12 E 2 /(l-y 1 y 2 ) , (3) 
D 12 = h 3 /12 E 1 y 2 /(l-y 1 y 2 ) = h 3 /l2 E 2 y 1 /(l-y 1 y 2 ) , (4) 
D 66 = h 3 /12 G , (5) 
where 


E^ = the effective modulus of elasticity (Young's 
modulus) in the ith direction, 
y^ = the effective Poisson's ratio in the ith 
direction , 

G = the shear modulus of elasticity (modulus of 
rigidity) , 

h = is the thickness of the plate. 

If the plate is simply supported on all four edges, the 
boundary conditions are zero displacement and zero bending 
moment along the edges. It can be shown that these two 
boundary conditions can be expressed as (see reference (6) ) 


5 



at = 0,d^ 


u = 0 , and 

(6) 

D n3 2 u/3cf 2 + D 12 3 2 u/3a 2 = 0 , 

(7) 


and at a 2 = 0,d 2 

u = 0 , and (®) 

D 12 3 2 u/3a^ + D 22 3 2 u/3a 2 = 0 . (9) 

The general solution of the homogeneous partial differential 
equation is given by a series expansion in the natural modes 
of the structure: 

u(a 1 ,a 2 ,t) = 

00 00 j (u) t+0 ) 

1 I n mn sin(mna 1 /d 1 )sin(mra 2 /d 2 )e mn mn . (10) 

m = i n= i 

where 

n = the influence coefficients of the modes in 
mn 

units of displacement, 

a) = the radian natural frequencies of the modes, 
mn 

and 

0 = the phase angles associated with the modes, 

mn 
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It is easily shown that this solution satisfies all of 
the simple support boundary conditions. Substituting this 
solution into the governing partial differential equation>of 
motion it can easily be shown that* the radian natural freq- 
uencies are given by the relation 


2 4, , 

“urn * " /ph 


D^Cm/d^' 


+ 2(D 12 +2D 66 )(m/d 1 ) 2 (n/d 2 ) 2 + D 22 (n/d 2 ) 4 


( 11 ) 


3.2 FREE VIBRATIONAL RESPONSE OF THE DAMPED CASE 

Including an equivalent viscous damping term, the 
governing equation becomes 


DjLl3 4 u/3aJ + 2(D 12 +2D 66 )3 4 u/3a 2 3a 2 + D 22 3 4 u/3a 4 

+Xu + phu = 0 . (12) 


By assuming that the vibrational modes of the plate are 
uncoupled, the modal damping coefficient can be defined as 

5 m „ = X/(2phii) mn ) . (13) 


Technically, the flexural modeshapes of a rectangular plate 
are uncoupled only when the damping is linearly related to 
the mass and stiffness properties of the plate in a special 
fashion. Thus, by making this assumption, our mathematical 
relations for the orthotropic plate become approximations. 
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The assumption that the modes are uncoupled can be used 
without serious error when the damping is a second order 
effect. Furthermore, the assumption is aided by the fact 
that the modal density for* flexural modes of a rectangular 
plate is nearly constant at low frequency. This means that 
there are no frequency regions of high modal density and 
therefore little opportunity for the exchange of vibrational 
energy among modes . 

Substituting into the governing equation, one obtains 

D n3 4 u/3aJ + 2(D 12 + 2D 66 )3 4 u/3a^3ci!2 + D 22 3 4 u/3a 4 

+ ph2E u) u + phu = 0 . (14) 

mn mn ' 

Assuming a solution of the form 
u(a 1 ,a 2 ,t) = 

oo oo yt+j0 

2 I n mn sin(m7ra 1 /d 1 )sin(mra 2 /d 2 ) e mn , (15) 

m= l n= i 

and substituting back into the governing equation, it can 
shown that the general solution to the damped free vibration 
problem is 
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u(a 1 ,a 2 ,t) = 


00 00 

5 I n ran sin(inira 1 /d 1 )sin(mra 2 /d 2 ) 
m= i n= i 


-Z 


0) t 
mn mn 


j(ft t+6 
J v mn mn 



( 16 ) 


where the damped radian natural frequencies are defined by 

o 1 1/2 


ft = to 
mn mn 


1-Z 


mn 


(17) 


Damping coefficients of the form given by equation (13) 
do not always accurately model the physical behavior exhi- 
bited by a plate. At higher frequencies the damping coeffi- 
cients of 2024 aluminum are relatively independent of fre- 
quency and are in the range of £ = 0.01. A more moderate- 
ly damped structure might have damping coefficients in the 
range 5 = 0.04. Because of the large energy loss due to 
sound radiation at low frequency, however, experimentally 
determined damping coefficients for the first few resonance 
frequencies will appear to be much, much larger than this. 
(On the order of £ = 0.1.) These inflated damping coef- 
ficients at low resonance frequencies are the result of 
neglecting the significant fluid loading effect in the model 
of the forcing function. For these reasons, two different 
damping models were investigated in this study given by the 
equations 
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( 18 ) 


^mn 


.01 

^mn 

• 06( “ll / %n ) + 

.04 


(19) 


where the hyberbolic term in each equation corrects for the 
fluid loading effect at the first few resonance frequencies 
and the constant term in each equation represents the inher- 
ent damping of the structure. Equation (18) shall be sub- 
sequently referred to as the small damping model while equa- 
tion (19) shall be referred to as the moderate damping 
model . 


3.3 FORCED VIBRATIONAL RESPONSE 

The general forced vibrational response of a simply 
supported, rectangular orthotropic plate is governed by the 

non- homogeneous partial differential equation 

i 


D lia 4 u/3oi! + 2(D 12 + 2D 66 )3 4 u/aa^9oi2 + 9 22 3 4 u/3a 4 
+ Xu + phu = q(a 1? a 2 ,t) , (20) 


where, for the purposes of this study, the forcing function 
q(a^,a 2 ,t) will be defined by the equation 

q( a ! , ° 2 ’ t) = 



6 ( Qf — c i ) 6 ( a 2“ c 2 ) + f a 



jwt 

e 


( 21 ) 
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The various parameters in equation (21) are 



the magnitude of the structureborne input 
(units of force), i 

the relative phase of the structureborne input, 
the magnitude of the airborne input (units of 
force/area), 

the relative phase of the airborne input, 

the location of the point vibrational input 

in the coordinate direction, i 

the location 'of the point vibrational input 

in the coordinate direction, 

the radian frequency of the simple harmonic 

forcing function, 


and the symbol 6 denotes the dirac delta function. Thus, 
the acoustic (airborne) input is chosen to be a normally 
incident, spatially uniform, simple harmonic forcing 
function and the vibrational (structureborne) input is 
chosen to be a simple harmonic, point vibrational load. 

The reasons for choosing a normally incident airborne 
input are twofold. First of all, the propeller noise which 
impinges on the sidewall of an aircraft defines an oblique 
incidence problem. The spatial pressure distribution of 
this incident sound is slowly varying due to the extremely 
long wavelengths of the low frequency sound. Therefore, 
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over small regions of the aircraft sidewalls, the spatial 
distribution of the pressure is nearly uniform. Secondly, 
it is comparatively much easier to construct an experimental 
apparatus that approximates the normal incidence condition. 

The point vibrational input was chosen to model the 
structureborne input primarily because it emulates the type 
of forcing function obtained experimentally using a shaker. 
It might be argued that a line load is a more realistic 
structureborne model since aircraft panels have ring frames 
and stringers attached to them. Because of the wide varia- 
tion in design of airframes and in wing attachment, however, 
it is not clear that a line load model would be any more 
realistic. Furthermore, for the purposes of understanding 
the interactive effects of the airborne and structureborne 
components, the point load model will, in principle, work 
just as well as a more complicated model. 

The forcing function of equation (21) can be expanded 
into a series representation using the eigenfunctions of the 
plate as follows: 


q(a 1 ,a 2 ,t) = 

r OO 00 

5 5 q mn sin ( miTa 1 / d 1 ) sin ( niT a2 //d 2^ 

- m = x n= i 


j(jdt 

e 


( 22 ) 


Utilizing the orthogonality principle it can be shown that 

the coefficients q m are given by 

mn 3 
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Utilizing the theory of complex variables, the term in 
brackets on the left side of equation (25) is given by 
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c ’ h “mn +ph25 mn ( 3“ ) -f >h “ 2 


phco 


mn 




1/2 


•^mn 
e , 


where the phase angle y mn is defined to be 


Y = tan 
1 mn 


-1 




Dividing equation (25) by equation (26) yields the 
participation factors 


-je. 


n e 

'mn 


mn 


s e + a e 

mn mn 


^mn 

e , 


where the modal influence coefficients are seen to 


s mn = 4 f s /( d i d 2 ph(JJ mn^ sin(mirc 1 /d 1 )sin(nirc 2 /d 2 ) 

1/2 


/ 


d-“ 2 /“L )2 + ' 2 5mn“ /u n, n )2 


and 


a mn - 4f a /tmn» 2 Ph M 2 n ) 


(1-cos (mir ) ) (1-cos (rnr) ) 


/ 


d-“ 2 /“mn> 2 + ' 2 5m„“/“ m r. )2 


1/2 


(26) 

(27) 
modal 

(28) 
be 

(29) 

(30) 
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8 U-'-vl n 


Thus, the steady state solution for the plate transverse 
displacement is 


CO CO 


u(a 1 ,a 2 ,t) = j j 


m= 1 n= 1 l 


s mn e + a mn e 


sin(miTa 1 /d 1 )sin(mra 2 /d 2 ) e 


jCuit-Y mn ) 


( 31 ) 


Simple differentiation with respect to time leads to an 
equation for the vibrational velocity of the plate given by: 

v(a 1 ,a 2 ,t) = 


I 


... i n= i l 


JO) 


-j<t>s -J4> a 

s mn e a mn e 


sin (mira 1 /d 1 ) sin (nira 2 /d 2 ) e 


j (ut- y ) i 


mn‘ 


( 32 ) 
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3.4 POWER DISSIPATED BY THE PLATE 

It can be shown that the time averaged power dissipated 
by the plate (see references [5] or [6]) is given by 


<W( t ) > . 


?2 ?! 


0 0 


1/2 Re{v*(r, t)v(r , t) } do^do^ 


d 2 ?1 


= X 


<v 2 (r,t]> t da^da 2 = XS<v 2 (r , t ) > r 


(33) 


0 0 


After lengthy complex algebra and straightforward 
integration, substitution of equation (32) into equation 
(33) leads to the result 


<w(t)> t = iu 2 /8 • 


I 5 

— 1 TV- 


xs 


ran 


mn 


2s a cos(<i> -6 ) 
mn mn VY s Y a' 


(34) 


The last term in equation (34) accounts for the interaction 
between the structureborne and airborne inputs in terms of 
the panel dynamics. Thus, equation (34) shows that the 
dynamic response of the panel to the combined airborne and 
structureborne inputs is, in general, not equal to the sum 
of the responses to the airborne and structureborne inputs 
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individually. Similar cross term components arise in the 
derivations of the power input to the plate. (See reference 
15].) One of the important features of equation (34) is 
that although it contains cross terms due to the combined 
inputs, it contains no cross terms between different modes 
of the structure. This absence of cross terms between modes 
is a direct consequence of the orthogonality principle 
between modes of the structure. Equation (34) indicates 
that when the structureborne and the airborne inputs are 
precisely 90 degrees out of phase, the inputs are 
uncorrelated and the cross term is zero. Thus, when the two 
inputs are uncorrelated, the power dissipated by the panel 
due to the combined inputs is exactly equal to the sum of 
the powers dissipated due to the structureborne and airborne 
inputs acting individually. 

4. ANALYSIS OF THE ACOUSTIC RESPONSE 
Once the surface velocity distribution of the plate has 
been found (equation (32)), the sound generated by the plate 
can be calculated. The classical approach for calculating 
the sound power radiated from a vibrating plate, which 
utilizes the principle of superposition of simple sources, 
was used in this study. This theory models each incremental 
area of the vibrating plate as a point monopole source near 
an infinitely rigid reflecting surface. The mathematical 
details of this classical theory are presented here for the 
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particular case of combined structureborne and airborne 
inputs. 

4.1 DERIVATION OF THE EQUATION FOR THE PRESSURE 

Utilizing Huygens principle of superposition of simple 
sources the sound generation of each incremental area, dS, 
of the plate is modeled as a monopole lying on a rigid 
planar reflecting surface. The incremental pressure in this 
case is given by (see reference [7]) 

dp ( r ' , t ) = 

- jkr ' 

jo)P o /(2irr ' ) e v(8 1 ,B 2 >t) dB 1 d8 2 • (35) 

In the 8^ system (see figure 1), the coordinates of the 
point of observation are (8^»3 2> 8 2 ) = ^ ^ l> r 2 ,r 3^ , an< ^ 
the coordinates of the incremental plate area are 
( B 1 * 8 2 * b 3 ) = (B 1 ,B 2 »0)- Thus, the 

i 

distance from a point in space to the center of the plate is 
given by 

r 2 = r 2 + r 2 + r 2 , (36) 

and the distance from the same point in space to incremental 
area on the plate which is radiating sound is 
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( 37 ) 


r> 2 = (r.-B,) 2 + (r 2 -B 2 ) 2 + r 2 

= r - 2r 1 6 1 - 2r 2 ^2 + + ^2 » 

where 


r^ = r sin(0)cos (<{>) , 

(38) 

r 2 = r sin(0)sin((f>) . 

(39) 

Substituting (38) and (39) into (37) 


r' = r 1-2 sin (©)cos (<f>) ( B-^/r ) - 2 

sin(0)sin(<})) (B 2 /r) 

2 2 1 1/2 
+ (8 L /r) 2 + (B 2 /r) 2 

(40) 

Using the far field approximations 


r >> d^ and r >> d 2 

(41) 


r' = 

r l-sin(0)cos(<j>) (B 1 /r)-sin(0)sin(<()) (B 2 /r) . (42) 

Substituting into equation (35), the expression for the 
incremental pressure caused by the incremental plate area is 

dp(r,t) = jwp o /(2Trr) 

- jkr( 1-sin (0) cos (<f>) ( B-j^/r )-sin (6) sin (<|>) (B 2 /r) ) 

• e 

• v(8 1 ,B 2 ,t) dB 1 dB 2 . (43) 
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Integrating over the area of the plate, the pressure at 

point in the far acoustic field is given by 

+ - jkr 

p(r,t) = ju)p o /(2irr) e 


d 2 /2 d x /2 


-d 2 /2 -d 1 /2 


jk sin (0) ( 3^cos (<j)) + B 2 sin (<{>) ] 


vCB-l, B lf t) 


dB 1 dB 1 


Making the following change of variables (see figure 1) 


Bf = Of - d 1 /2, B 2 = o<2 “ d 2 //2 


dBjL = da 1? 


dB 2 = da 2 , 


equation (44) becomes 


P ( r , t ) = jtop o /(27rr) 


-jk(r + d^/2 sin (0)cos (<f>) + d 2 /2 sin(0)sin(<}>) ) 


d d 

i* ^ r r- 


0 0 


v(a 1 , a 2 , t) 


jk sin(0)(a 1 cos(<j>) + a 2 sin(<}))) 
e 


da^da 2 . 


(44) 

(45) 

(46) 


(47) 


20 



Returning to equation (32) and substituting equation (28) 
into (32), the velocity of the plate can be shown to be 

i 

jwt 

v(a 1 ,a 2 ,t) = j(j e 

co co _ je 

• I 1 . ^mn e mn s i n ( m7Ta i/ d i ) s i n ( n7ra 2 //d 2^ ’ (48) 

m = i n= i 

Substituting equation (48) into equation (47) 
p(r,t) = -w 2 p o /(27rr) 


j(wt-k(r + d^/2 sin(©)cos (<}>) + d 2 /2 sin (©) sin (<{>)) ) 
• e ^ 


I I 

m = i n = i 


-j© 


^mn e 


mn 


1 

jk sin (Q)cos (<})) 

e sin (imra^/d^ ) da^ 

0 


d„ 

■ z jk sin (0) sin (<j>) a 2 

e sin(nira 2 /d 2 ) da 

0 


(49) 
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4.2 EVALUATION OF THE INTEGRALS 

The integrals in equation (49) can be written as 


, d i 


m 


jk sin(6)cos (<}>) a-^ 


sinCmira^/d^) da^ 




cos(k sin (0)cos (<j>) a-^) sin (mua^/d^ ) da^ 


sin(k sin (e)cos (<}>) a^) sin (mira-^/d^ ) da^ . (50) 


and 


?2 


jk sin(©)sin(<j>) 


sin (niTa 2 /d 2 ) da 2 


?2 


cos(k sin (© ) sin ( <f>) a 2 ) sin(mta 2 /d 2 ) da 2 


sin(k sin(0)sin(<{)) a 2 ) sin(ntra 2 /d 2 ) da 2 . (51) 


With the integrals represented in this form, they can be 
evaluated by making a trigonometric substitution followed by 
straightforward integration, or more simply by looking them 
up in tables of integrals. 
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4.3 CALCULATION OF THE SOUND POWER RADIATED 


With the integrals of equation (49) evaluated as dis- 
cussed in the previous section, the equation for pressure 
can be written as: 

p(r,t) = - 0 ) 2 p o /( 2 irr) 

j(wt-k(r + d 1 /2 sin(0)cos(<j)) + d 2 /2 sin (0) sin (<j>) ) ) 


OO 00 

• I I 

m= i n= i 


ran 




(52) 


From this last equation, the acoustic pressure at an arbi- 
trary point in space can be calculated where it has been 
assumed that the point of observation is in the far acoustic 
field, i.e. r >> d^ and r >> . The time averaged 

magnitude of the acoustic intensity vector in the radially 
outward direction at a far-field point in space for a simple 
harmonic, monopole source is given by 

<I r ,(r',t)> t = |p(r' , t) | 2 /(2p Q c o ) . (53) 
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This relationship also holds for the simply supported 
rectangular plate when r >> d^ and r >> d^ • Making the 
appropriate substitutions, the far field time averaged 
magnitude of the acoustic intensity vector in the radially 
outward direction at a point in space for the plate 
vibrating in simple harmonic motion is given by 


<I r (r,t)> fc = p (r, t)p(r , t)/(2p Q c o ) 


(p o u) 2 /(2iTr ) ) 2 /(2p o c o ) 


OO 00 

L] S J 

L k= l 1=1 


+ je 


kl * * 

n,,! e z v z 


kl 


k “1 


OO OO 

5 I 

m= i n= i 


^mn 

^mn e Z m Z n 


(54) 


Rearranging the order of summation 

.2 / r ^ \ 


oo oo oo oo 


<1 


r ( r , t ) > fc = (p Q m /(2irr)) /(2 p o c q ) 5511 

k = i 1= i m= i n = i 


+ -^ e kl ^ 6 mn * * 

n k i e 1 %n e Z,_ Z, Z_ Z 


Ji U 1 U L. I 

k 1 m n 


(55) 


The time averaged sound power is given by 


<TT( t ) > , 


<I*n> t dS . 


(56) 
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In the present case, the surface to be integrated over 
hemisphere in the far acoustic field, and the normal u 
vector is in the radially outward direction. Thus 


<TT( t ) > . 


2ir ir/2 


0 0 


<I r (t)> t r 2 sin(6) d0d<j> . 


The time averaged sound power radiated by. the plate is 
therefore given by 

<7T( t ) > , = 


2 tt ir/2 


O' 0 


OO 00 00 00 


(p o w 2 /(2it)) 2 sin(e)/(2p Q c 0 ) 1111 

k= i l=i m= i n = i 


+j6 kl - j0 mn * * 

n, , e g e Z, Z, Z Z 

'kl 'mn k 1 m n 


dedcj> . 


From equation (28) it can be shown that 


+ i®kl -3®mn ^kl'W 

"kl e "mn e - e 


s, , s + a, , a 
kl mn kl 


MW . 

+ 3 kl a mn e + a kl s mn e 


is a 
it 

(57) 


(58) 


mn 

(59) 
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Substituting equation (59) into equation (58) the result is 


<TT( t) > t = 


2 it tt/2 

f f r oo oo oo oo 

(p o w 2 /(2tt) ) 2 sin(e)/(2p o c o ) I ^ I ^ 

J J L k= i 1= i m= i n= i 

0 0 


kl mn 


+ a 


kl mn 


+ s 




kl a mn 


+ a 




kl mn 


^ Y kl Y mn^ 

• e 



d0dcj> . 


(60) 


Since the sound power is a real quantity, the series in 
equation (60) must be real. Using simple complex algebra, 
the equation for sound power can be written as: 


26 



<TT(t) >, 


2 it tt/2 

r /• 


o o 
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d0d<j) 


(61) 


The last term in equation (61) accounts for the interaction 
between the structureborne and airborne inputs in terms of 
the sound power radiated. Thus, equation (61) similarly 
shows that the sound power generated by the panel due to the 
combined airborne and structureborne inputs is, in general, 
not equal to the to the sum of the individual sound powers 
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radiated by the panel due to the airborne and structureborne 
inputs acting individually. Unlike equation (34) for the 
panel dynamics, however, equation (61) contains cross terms 
between the different modes of the structure. The addi- 
tional cross terms between different flexural modeshapes 
arise because the spatial average of the acoustic intensity 
used in calculating the sound power takes place over a hemi- 
spherical surface in the far acoustic field, and not over 
the surface of the plate (the area of spatial averaging for 
the analysis of the dynamic response of the plate). Thus, 
the orthogonality principle does not apply when calculating 
the sound power radiated by the plate. Furthermore, 
equation (61) indicates that even when the structureborne 
and airborne inputs are precisely 90 degrees out of phase 
(uncorrelated), there are still many non-zero cross term 
components between the two inputs which contribute to the 
total sound power radiated due to the existence of the cross 
terms between different flexural modes of the structure. 
Therefore, in general, the total sound power radiated due to 
the simultaneously combined structureborne and airborne 
inputs is never exactly equal to the sum of the sound powers 
radiated due to the structureborne and airborne inputs 
acting individually. 
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5. RESULTS 


This section contains selected results of a numerical 
implementation of the analytical model. Source code list- 
ings of the computer programs used in the analytical studies 
are given in reference 1 5.] . The analytical results shown 
here are intended to simulate the noise radiation of a . 8 mm 
(0.032 inch) thick plate constructed of AA 2024 aluminum 
(surface density of 2.22 kg/m ) due to both the individual 
and combined acoustic and vibrational inputs. The physical 
dimensions of the plates were chosen to be 0.406 m x 0.241 m 
(16 in x 9.5 in) since an extensive study of the transmis- 
sion loss properties of plates of this size had already been 
completed by NASA researchers. Unless explicitly stated 
otherwise, the results are for the case of a normally 
incident, spatially uniform, 1 Pa peak acoustic load and a 
0.01 N peak point vibrational load located at the coordi- 
nates of a^= 0.060 m and a 2 = 0.135 m. Also, unless 
stated otherwise, the small damping model (equation (18)) 
was used in the simulations. In all cases the forcing 
functions were applied uniformly over the frequency range at 
2 Hz intervals. 

The results that follow are divided into three sec- 
tions. In the first section, the effects of the relative 
magnitudes and phase of the inputs are investigated. In 
second section, the effects produced by changing the shaker 
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location (altering the path of the structureborne input) are 
examined. Finally, the third section examines the effects 
of adding damping treatment to the aluminum plate. 

5.1 EFFECTS OF RELATIVE MAGNITUDE AND PHASE 

5.1.1 Analytical Results for a Dominant Airborne Source 

Figure (2) shows the sound power levels produced by the 

individual airborne and structureborne inputs. This figure 
shows that, in this case, the airborne source is dominant. 

The sound powers produced by combining the airborne and 
structureborne inputs with phase angle differences of 90 
degrees, 0 degrees, and 180 degrees are shown in figures 
(3), (4), and (5) respectively. Figure (3) shows that the 
sum of the results of the individual inputs is roughly 
equivalent to results obtained by combining the inputs at 90 
degrees over most of the frequency range. Figure (3) also 
shows that, at least in some isolated frequency regions, 
,ther^ can be a considerable difference between the sum of 
the individual sound power components and the combined sound 
power of the uncorrelated inputs. Examination of equation 
(61) leads one to the conclusion that the differences in the 
two curves of figure (3) can only be caused by cross terms 
between the inputs that represent cross coupling between 
different modes of the structure. Figures (4) and (5) show 
that even larger deviations between sound potoers produced by 
the sum of the individual inputs and the combined inputs 
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occur when the inputs are combined at 0 or 180 degrees. 

These results underscore the importance of the cross terms 
contained in equation (61) by showing that their contri- 
bution to the overall sound power radiated can be as large 
as the dominant term. Furthermore, since the sign of the 
cross terms can be positive or negative, the overall noise 
generated can vary over a very large range depending on the 
relative phase of the inputs. 

5.1.2 Analytical Results for a Dominant Structureborne 
Source 

Figure (6) shows the computed sound power levels pro- 
duced by the aluminum plate due to a 0.2 Pa peak acoustic 
load and an independent 0.05 N peak vibrational load. This 
figure shows that, in this case, the structureborne source 
is dominant. 

The sound powers produced by combining these two inputs 
with phase angles of 90 degrees, 0 degrees, and 180 degrees 
are shown in figures (7), (8), and (9) respectively. Part-, 
icularly good agreement between the sum of the individual 
inputs and the case of uncorrelated, combined inputs is seen 
to occur in figure (7). This is an indication that the 
cross coupling between different modes of the structure play 
a relatively minor role in this case. Figures (8) and (9) 
again shows that large deviations between the cases of the 
sum of the individual inputs and the combined inputs occur 
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when the inputs are combined at 0 or 180 degrees. Thus, the 
cross terms between the inputs for like modes can still be 
responsible for a significant portion of the overall sound 
radiation in this case. 

5.2 EFFECTS OF ALTERING THE INPUT PATHS (SHAKER LOCATION) 
5.2.1 Analytical Results for a Dominant Airborne Source 

Figures (10) through (13) show the results that were 
obtained from the analytical programming for the case of a 
normally incident uniform acoustic load and a point vibra- 
tional load located at the new coordinates of a^= 0.121 
m and 0.203 m. These coordinates place the vibra- 

tional point load at the center of the plate, thus driving 
the same modes as the acoustic input (i.e. the odd modes of 
the plate). 

Figure (10) shows the sound power levels produced by 
the aluminum plate due to a 1 Pa peak acoustic load and an 
independent 0.01 N peak vibrational load. This figure shows 
that the airborne source is dominant or at least of equal 
influence over the entire frequency range for these inputs. 

The sound powers produced by combining the two inputs 
are shown in figures (11) through (13). These figures show 
that the overall noise radiation due to combined inputs is 
extremely sensitive to changes in the relative phases 
between the inputs. This result indicates that the cross 
terms have an even larger influence on the overall noise 
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radiation if and when the airborne and structureborne inputs 
drive the same modes of the structure. 

5.3 EFFECTS OF ADDED DAMPING 

5.3.1 Analytical Results for a Dominant Airborne Source 

Figures (14) through (17) show the results that were 
obtained from the analytical program for baseline conditions 
(shaker located at 0.060 m and « 2 = 0.135 m) on 

the aluminum plate using the more moderate level of damping 
included in the model (equation (12)). 

Figure (14) shows the sound power levels produced by 
the aluminum plate due to a 1 Pa peak acoustic load and an 
independent 0.01 N peak vibrational load. This figure shows 
that the airborne source is highly dominant over the entire 
frequency range. 

The sound powers produced by combining these two inputs 
are shown in figures (15) through (17). The sound power 
curves of figure (15) suggest that cross coupling between 
different modes occurs in the same frequency regions as 
figure (3) (the case of small damping). Figures (16) and 
(17) show that the overall noise radiation due to combined 
inputs is again very sensitive to changes in the relative 
phases between the inputs. The figures show that the sound 
power level can vary as much as 5 dB in isolated frequency 
regions due to the influence of the cross terms. The 
increased smoothness of the sound power curves, due to the 
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added damping, makes the effects of the cross terms parti- 
cularly evident. 

6. DISCUSSION AND CONCLUSIONS 
To the author's knowledge, this is the first analytical 
study of the problem of the sound radiation of aircraft type 
materials due to fully coherent combined airborne and struc- 
tureborne inputs in the low frequency regime. The single 
most significant finding in the results of this study was 
the discovery of the relative importance of the cross term 
sound power components. The importance of this finding can 
not be over emphasized in light of the fact that all of the 
previous analytical and experimental studies performed in 
the last 7 or 8 years have neglected the interaction between 
the airborne and structureborne inputs in the low frequency 
regime. In a recent review paper (see reference [1]) it was 
found that 93 publications appeared in the open literature 
between 1978 and 1984 which devoted their attention strictly 
to the airborne noise transmission through the sidewalls of 
aircraft. A smaller number of studies (see references [2]- 
[4]), addressing the problem of structureborne noise paths 
in aircraft, similarly restricted the scope of their anal- 
ysis. As a result of this study, some of the conclusions of 
these earlier studies regarding the effectiveness of certain 
noise control treatments may now be in doubt. 
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Evidence to support this claim can be found in the 
results such as those found in figures (14) through (17). 
Computations on the analytical results in these figures (for 
the damped aluminum panel) show that while the airborne 
component contributes 82.0 dB overall and the structureborne 
component contributes 68.6 dB overall, the levels for the 
combined sound power ranged anywhere from 80.7 dB to 83.3 dB 
overall. Thus, the overall level of the combined sound 
power varied over a 2.6 dB range depending on the phase 
relationship between the airborne and structureborne inputs. 
Simple calculations show that the sum of individual airborne 
and structureborne components is 82.2 dB or, in other words, 
one would expect a 0.2 dB increase in the overall sound 
power level due to the addition of the structureborne noise 
to the dominant airborne noise if the sound powers were 
additive. The implication here is that the structureborne 
component does not have to radiate a significant amount of 
noise on its own in order to significantly change the level 
of the combined sound. In fact, additional results pre- 
sented in reference 1 5 ] indicate that the structureborne 
component needs only to change the dynamic response of the 
structure in order to also significantly influence the over- 
all noise radiation. 

The results of the study alSo indicate that the inter- 
active effects are due primarily to cross terms between like 
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modes of the airborne and structureborne inputs. While the 
results for uncorrelated inputs showed that the effects of 
cross coupling between different modes of the structure 
could be significant in some frequency ranges, the largest 
interactive effects were found when the airborne and struc- 
tureborne components were correlated. It should also be 
noted that, since the cross coupling of different modes 
enters the problem through the damping model (see equation 
(61)), the contribution to the sound power due to cross 
coupling between different modes of the structure may have 
been a spurious effect that was introduced through the use 
of equations (18) and (19). Thus, the effects on the com- 
bined sound powers due to cross coupling between different 
modes may or may not be observed experimentally. 
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Figure 1 - Geometry of the simply supported ■ plate . 
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Figure 2 - Predicted sound powers of the individual components for a dominant 
airborne source and the shaker located near the corner. 
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Predicted sound powers for the case of the inputs combined at 0 degrees 
with a dominant airborne source and the shaker located near the corner. 
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Figure 6 - Predicted sound powers of the individual components for a dominant 
structureborne source and the shaker located near the corner. 
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Figure 10 - Predicted sound powers of the individual components for a dominant 
airborne source and the shaker located at the center. 
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Figure 12 - Predicted sound powers for the case of the inputs combined at 0 degrees 
with a dominant airborne source and the shaker located at the center. 
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Figure 14 - Predicted sound powers of the individual components for a dominant 

airborne source, the shaker located near the corner, and a moderate level 
of damping added to the panel. 
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Figure 15 - Predicted sound powers for the case of the inputs combined at 90 degrees 
with a dominant airborne source, the shaker located near the corner, and 
a moderate level of damping added to the panel. 
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Figure 16 - Predicted sound powers for the case of the inputs combined at 0 degrees 
with a dominant airborne source, the shaker located near the corner, and 
a moderate level of damping added to the panel. 
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Figure 17 - Predicted sound powers for the case of the inputs combined at 180 degrees 
with a dominant airborne source, the shaker located near the corner, and 
a moderate level of damping added to the panel. 
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